Effect of nuclear interactions of neutral kaons on CP asymmetry measurements 



O 

(N 

> 
O 

(N 



X 

CD ' 
i ' 

43 . 



(N 
> 
oo 
m 
o\ 

o 
o 



X 



B. R. Ko, 1 E. Won/'B B. Golob, 2 ' 3 and P. Pakhlov 4 

1 Korea University, Seoul 
^Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana 
3 J. Stefan Institute, Ljubljana 
'''Institute for Theoretical and Experimental Physics, Moscow 

We examine the effect of the difference in nuclear interactions of K° and K° mesons on the 
measurement of CP asymmetry for experiments at e + e~ colliders - charm and B-meson factories. 
We find that this effect on CP asymmetry can be as large as 0.3%, and therefore sufficiently 
significant in interpreting measurements of CP asymmetry when neutral kaons are present in the 
final state. 
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Modern high-statistics B factories discovered the joint 
violation of charge-conjugation and parity {CP) in B- 
meson decay modes. In some B° decays Q, large CP 
violation induced by B° — B° mixing is observed to be 
consistent with the predictions of the standard model 
(SM) and the Kobayashi-Maskawa ansatz Smaller, 
direct CP violations, attributed to the interference of 
different amplitudes, but without mixing have also been 
reported 0, SM predictions for the direct CP vi- 
olation in many charmed-meson decays are typically of 
However, the present accuracy of measure- 
ments of CP asymmetry in D meson decays is close 
to their SM expectations. For example, in the decay 
D + KgTr + [6|], the statistical sensitivity on the mea- 
sured CP asymmetry (of ~0.2%) [7| is slightly smaller 
than the effect expected in the SM of (0.332±0.006)% 
from K° — K° mixing Q. Experiments at future high- 
luminosity B factories and at the LHC are likely to reach 
the sensitivity needed to observe CP violation in some 
D decay modes. 

The measured asymmetries of B or D mesons for de- 
cays which has Kg in their final states, can be mimicked 
(or diluted) by differences between K° and K° interac- 
tions with detector material. The probability of an in- 
elastic interaction of a neutral kaon in the detector de- 
pends on the strangeness of the kaon at any point along 
its path, which is due to oscillations in kaon strangeness 
and different nuclear cross sections for K° and K° . Hence 
the total efficiency to observe a final state K s differs 
from that expected for either K° or K°. This effect is 
related to the coherent regeneration of neutral kaons Q . 
This kind of contribution may be non-negligible for pre- 
cise measurements of direct CP violation in B and D 
decays, and also important in the determination of fa 



in the B 



D°K^ 



(K s ir + ir )dK + transition 



distorted by the K° interaction. 

In this paper, we evaluate the effect of the differ- 
ence in nuclear interactions of neutral kaons on mea- 
surements of CP asymmetry performed at charm and B 
factories, or will be carried out at the near future high- 
luminosity B factories. Our study represents an exten- 
sion and more detailed description of the method used 
to estimate the effect of K° /K° interactions in material 
in Ref. Q- We also note that the detector-simulation 
program Geant4 commonly used in high energy 



physics experiments, does not take into account the ef- 
fect considered in this paper, as the K° and K° are pro- 
jected onto the Kg or K\ components at their production 
point rather than at their points of 7T7t decay. The time- 
dependent K a — K° oscillations are thereby ignored in 
Geant4. A similar effect in D° — D° oscillations was 
found to be small in the mass and lifetime differences 
between D° and D° [e| • The aim of this paper is to ap- 
proximately estimate the magnitude of the effect due to 
the difference in K° and K° nuclear interactions under 
conditions of current and future experiments, and bring 
this issue to the simulation developers for possible in- 
clusion in programs such as Geant4. The method and 
result can serve as an estimate of systematic uncertainty 
for measurements neglecting the effect, or as a starting 
point for more refined calculations to be used in the fu- 
ture experiments in order to correct for the effect. 

Let us consider production of some meson V and its 
antimeson V in e + e~ collisions, each followed by its de- 
cay into states containing a neutral kaon, and observed 



through the K s 



V 
V 



K° s 



mode: 

X, 
X. 



and in a precise measurement of D° — D° mixing in the 
KgTr + ir~ final state, as the Dalitz distribution would be 



V can be a charmed or B meson. For certain charmed 
meson decays, there is a small contribution from doubly 
Cabibbo-suppressed decays that we ignore, in our main 
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calculation, but assign a systematic uncertainty for this 
assumption. The CP asymmetry in the V decays is de- 
fined as 
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CP 
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(1) 



where T denotes the partial decay width. We assume 
that the production point is surrounded by a cylindrical 
structure of material, typically used in a collider detector 
environment, such as a beam pipe and several thin layers 
of vertex detectors. 

To obtain the time development of neutral kaons in 
matter, we use the calculation carried out in Refs. 
14| . The time evolution of amplitudes in the K^ and Kg 
basis, as given in Ref. [14| . becomes 
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The quantities cei,(t) and as(^) are the amplitudes for 
finding states as a if£ and i-Tlj at some time respec- 
tively, and a° and ag are those states at t=0, where 
t refers to their proper times. The masses and ms, 
and decay widths Tl and T$ refer to and K s , respec- 
tively. The quantity m in A\ denotes the mass of the 
K° and K° . The volume density of material is Af= ^jf- , 
where p is the mass density. Na is Avogadro's number, 
and M is the mean molar mass. The quantities / and 
/ are the forward scattering amplitudes of K° and K°, 
respectively. The parameter r is called the regeneration 
parameter, defined as r = and its magnitude is gen- 

erally small, typically in the order of 10 -2 . Expanding 
a^(t) and as(t) up to the first order in r, we obtain 



a L (t) = + C(*)«sr, 

a s (t) = &(*)<*§ + C(*Kr, 
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where %L,s(t) = 
^q— *Az,i 
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From these relations, the 



amplitudes following the passage of several layers of de- 
tector material, can be obtained iteratively as follows: 

a h(tj) = - tj-i)ah(tj-i) 

<*s{tj) = €s(tj -tj-i)a s {tj-i) 

+ C(tj-fj-i)aL(*i-i)rj (4) 

where index j refers to the layer of material last pene- 
trated. This follows because neutral kaons pass through 
several layers of vacuum and material before they de- 
cay. The number of terms in Eq. §Q increases rapidly as 
the number of detector layers increases, and it is squared 
when |aL,sW| 2 are computed to obtain the probability. 
We evaluate all terms using the symbolic calculation pro- 
gram Mathematica Our dilution effect (A-p) can 
be extracted from the total asymmetry (At), which is 
incorporated in K° regeneration. Without CP violation 
in V — s- Kg + X decay itself, the At in the decay can be 
expressed as 

jR(t) dT v ^ K s+ x -JR(t) dT p ^ K s+ x 
J R(t) dT v ^ K °s+ x + J R(t) dT p ^ K °s+ x 
A v + A 
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where R(t) and R(t), the two-pion decay rates for initial 
K° and K°, respectively, can be expressed as Rs\as(t) + 
770^ (i)| 2 - Rs is the time-independent decay rate of the 
Kg eigenstate, and the ratio of amplitudes r\ — M. (K\ — > 
h+tt-)/M(K s -> tt+tt-). The first term in Eq. ©, 
Aq P: is the asymmetry due to K° — K a mixing which 
is not of primary interest in this paper, and thus can be 
subtracted. The third term, A[ n t is the asymmetry from 
interference between the CP violation in K° mixing and 
the material related amplitudes, and is expected to be of 
0(|r?7|) ps 10~ 5 . We estimate the third term numerically 
as fa 10 -6 , and therefore ignore it. Hence, At reduces 
to A v if the CP violation effect due to K° - K° mixing 
in At is removed, thereby setting the parameter 77=0. 
Approximating *He(A/)/3m(A/) = 1, and Am 
A-p can be expressed as 
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A v oc - 3m(A/) oc <t(K°N) - a(K°N), 



(6) 



where N refers to the atomic nucleon of the detector 
material. 

To compute results for Eq . (|6]) taking into account ef- 
fects of nuclear screening [16| , we adopt an empirical scal- 
ing law based on measurements in C, Al, Cu, Sn, and 
Pb for neutral kaon momenta (p ) between 20 and 140 
GeV/c 
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FIG. 1: The Acr(K~{d,p}) = a(K-{d,p}) - a{K + {d,pV) 
values as a function of kaon momentum, obtained from 
for the proton (dotted lines) and the deuteron data (solid) 
are shown in the top plot. The ratio of two cross section 
differences as a function of kaon momentum is shown in the 
bottom plot (solid circles), together with the fit using the 
error function (curve). 



where A is the atomic number and 0.758 accounts for 
nuclear screening. The scaling of A - 758 in Eq. ([7]) also 
describes Pb, Cu, and C data quite well^down to 5 
GeV/c The deuteron data in Ref. 

well with the prediction of Eq. ([7]) for A 
to 200 GeV/c. We extend the scaling down to lower 



18j also agree 
2 from 50 



momenta assuming isospin symmetry of nuclear inter- 



actions, a(K°ii)=a(K--p) and a(K°p)=a(K^ 



We 



approximate a(K + n)=a(K + p) to improve the estima- 
tion of A-p, and this assumption is consistent with mea- 
surements Q. (Symbols p and n correspond to the 
proton and neutron, respectively.) Using experimen- 
tal results for a{K~p) and a(K + p) from Ref. 0, we 
obtain Aa(K~ {d, p}), where d denotes deuteron, with 
Aa(K-d)=a(K-d)-<j(K+d) and Aa(K-p)=a(K-p)- 
a(K + p). Figure [1] shows Aa(K~ {d, p}) (top) and the ra- 
tio of the two, Aa(K~ d) / Ao(K~ p) (bottom), as a func- 
tion of the kaon momentum. We fit the ratio of Aer( K~d) 
to Aa{K~p) using an empirical function while keeping 
the nuclear screening term A 758 fixed. The value of 
X 2 /d.o.f is approximately 2, indicating our modeling of 
the ratio of Aa(K~d) to Aa(K~p) is not unreasonable, 
so that Eq. (UJ, obtained in the high- momentum range, 
can be scaled down to 1 GeV/c and below. Using the fit, 
Eq. (J7J is altered as follows 



Aa{K°N) = 



A°- 75S Aa(K-p) 



1 + I.252e- 1 - 841 P A ""( GeV / c ) 



mb, (8) 



where p K is the momentum of K . We use Eq. © 
in the numerical calculation of Eq. ([5]). The numera- 



tor in Eq. ((5J should extrapolate the screening effect to 
atoms in the detector material we use in Table U and 
the denominator reflects the low-momentum behavior of 
the difference in cross section between the proton and 
deuteron data. We compared our scaling method with 
the experimental data [19j, and found a good agreement. 



> < 

C5 , 




4 6 

p K ° s (GeV/c) 



^ 8000 L 



CD 

" 6000- 

f^ooo: 

I 2000- 
■ ■I 




08000 
d 



4 6 

p K ° s (GeV/c) 



jg 8000 - 
<D 

CD 

6000- 

p 

2-4000- 

<5 2000- 
LU 



O 8000 

o 

« 6000 
W 4000 
2000 



1 -0.5 


0.5 

cose 


: :B°- 

: :B°- 


>K° s n + 

> K° s n + n- 

> KgK + K" 




1 -0.5 


0.5 K 
cose 


: ;d°- 

: D°- 


>K°7I + 

> KgK + K" 




1 -0.5 


cose 



p s (GeV/c) 

FIG. 2: Kg momentum (left column) and angular distribu- 
tions (right column) for different decay modes. The upper two 
rows are for y / s=I0.58 GeV and the ,07=0.425 configuration 
and the bottom row is for -^5=3770 MeV. 

To obtain the expected four-vectors of K s mesons in 
the final state, we use Pythia ^20] and EvtGen 21 1 
Monte Carlo codes to simulate generation and decay of 
charmed and B mesons produced in e + e~ collisions. Two 
kinematic cases are considered reflecting two distinct ex- 
perimental environments: the first case is for a center- 
of-mass energy y / 5=f0. 58 GeV and a Lorentz boost fac- 
tor of /?7=0.425 (B factory), and the second case is for 
i/s=3770 MeV with no Lorentz boost (charm factory). 

The numerical values of Eq. © are calculated for 



K^+, 



B { > 



K%TT+- 



K%K+K~ 



K%K+K- 



7i- , -D u 
and B° 

produced in the two kinds of e + e~ collisions described 
above. The choice of the decay channels is arbitrary, but 
intended to show a broad range of momenta that depend 
on decay characteristics. The first four plots in Fig. [2] 
show the momentum and polar angle distributions of K s 
mesons in the laboratory frame for served final states at 
Vs=10.58 GeV and /37=0.425. The distributions in po- 
lar angle are seen to be very similar for K s from charmed 
and B meson decays, despite that the momentum distri- 
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TABLE I: Two beam pipe and detector configurations se- 
lected for the study described in the text, with S and R cor- 
responding to the thickness and radius of the given detector 
component. There are two configurations of layers given for 
Case II. 



Material 


Beam pipe 
Be 


Detector layers 
Si 


Case I 


(5=1 mm 

at R—1.5 cm 


5=300 urn. 

at i?=2.0, 4.35, 7.0, 8.8 cm 


Case II 


5—1 mm 

at R=1.0 cm 


<5=50 fim 

at R=1.4, 2.2 cm 

(5=300 /um 

at i?=3.8, 8.0, 11.5, 14.0 cm 



butions show large differences among the decay modes, 
which causes significant differences in the values of Ad- 
As for the material geometry, we choose two general 
detector options, summarized in Table HI that closely re- 
semble the existing or planned B-meson and charm fac- 
tories. The first option, denoted as "Case I" [22], |23j . 
reflects the current charm and i?-meson-factory experi- 
ments. The second option, denoted as "Case II" , reflects 
a proposed super B- factory experiment [24j. We apply 
typical geometrical acceptance criteria in calculating Ad 
for each case. 

We calculate Ad for Case I, with y/s— 10.58 GeV and 
,37=0.425, for the decay modes mentioned previously, 
and their resultant values are summarized in Table HU 
We find that Ad values are «10 -3 for all the above decay 
modes, and they are mainly affected by the beam pipe. 
We also plot the distributions of Ad as a function of mo- 
mentum and polar angle of Kg for Case I. The upper 
plots of Fig. [3] are the Ad distributions for D + — > Kgir + 
at ^5=10.58 GeV and /3 7 =0.425. The values of A v de- 
pend strongly on Kg momentum distributions as shown 
in Fig. [3] and are larger for smaller momenta. This can be 
understood from the fact that the cross section difference 
is larger at small momenta as shown in the upper plot 
of Fig. [TJ We apply typical experimental selection crite- 
ria of p*{D+) > 2.5 GeV/c and pt{tt + ) > 0.45 GeV/c 
in D + K° s ti + decay, where p*(D + ) and pr(i" + ) are 
the momenta of D + in the center-of-mass frame and the 
transverse momenta of tt + in the laboratory frame, re- 
spectively. We find practically no difference in Ad ap- 
plying these selection criteria. 

The major systematic uncertainty in this calculation 
is from the assumption £He(A/)/ 3m(A/)=l. We es- 
timate this effect using momentum-dependent values 
of 9\t(Af)/ 3m(A/), where ffte(Af) is obtained from 
the best known values in Ref. |25j (kaon momenta 
available up to 2.6 GeV/c). The results differ from 
«He(A/)/3m(A/)=l by 6% when limiting the K% mo- 
mentum range up to 2.6 GeV/c for charmed-meson de- 
cays. For B meson decays, the effect is found to be 10%. 
Because of limited information on CHe(A/), we assign 



systematic uncertainties of 10% and 20% for charmed 
and B meson decays, respectively, for the assumption of 
£Re(A/)/ 3m(A/) = l. The systematic effect from the as- 
sumption that Am ~ \ AT is found to be negligible. Sys- 
tematic effects due to uncertainties in modeling Eq. © 
are also found to be negligible. The systematic uncer- 
tainties from the measurements for a(K~p) and a(K + p) 
are 0.5% and 0.9%, respectively. Systematic uncertain- 
ties due to the statistical uncertainties on a{K~p) and 
a(K + p) are estimated from Monte Carlo, and found to 
be negligible. Other sources include uncertainties on Am, 
and lifetimes of K\ and Kg, and are also negligible. 
There is a contribution from doubly Cabibbo-suppressed 
decays of charmed mesons that is neglected in the com- 
putation of Ad- According to Ref. |2(|, we assign a 10% 
systematic uncertainty to the final states with a contri- 
bution from doubly Cabibbo-suppressed decays. 

In the study of the same decay channels of charmed 
mesons for the center-of-mass energy in the region of 
V>(3770), we introduce no Lorentz boost for the detec- 
tor geometry described by Case I. This checks the effect 
of different kinematics of Kg by comparing the results 
with those for ^=10.58 GeV and ^7=0.425. The bot- 
tom two plots in Fig. [2] show the momentum and po- 
lar angle distributions of Kg mesons for ^=3770 MeV 
with no Lorentz boost, showing lower Kg momentum dis- 
tributions relative to those from the configuration with 
y / s=10.58 GeV and /3 7 =0.425. Larger Ad values are con- 
sequently expected, which is consistent with the calcula- 
tions shown in the third column of Table Ql] The bottom 
of Fig. [3] shows the distributions of Ad as a function of 
momentum and polar angle of Kg for D° — > KgTr + ir~ at 
Vs=3770 MeV and /3 7 =0. We find that the A v values 
are in general larger than given in the second column of 
Table HI1 Again, this reflects the Kg momentum distri- 
bution shown in the bottom plot in Fig. [5J which peaks 
in the phase space region with the largest Aa(K° N). 
Here, the systematic uncertainty from the assumption 
that 9te(A/)/3m(A/) = 1 is found to be 30%, and other 
sources are negligible. 

As a final benchmark, we also evaluate Ad for the Case 
II configuration with ^=10.58 GeV and /3 7 =0.425. The 
results of estimating of Ad are listed in the last column 
of Table [III This configuration checks the effect of differ- 
ent geometry for detector material by comparing results 
from Case I with the same kinematics. We find that the 
contribution of the first two thin layers of Si sensors is 
negligible. Furthermore, the contribution of the outer Si 
sensors is also smaller than that of the Si sensors in Case 
I as their distances from the production point of neutral 
kaons are longer. This results in smaller dilution than 
for Case I. Systematic sources and effects are similar to 
those of due to Case I. 

As shown above, the dilution effect in the calculation of 
Ad is most sensitive to the momentum of Kg, and mainly 
due to the beam pipe contribution. Hence, the dilution 



5 



TABLE II: Numerical estimation of Ad for three configurations. The values in parentheses are only for the beam pipe element. 







Configurations 






Case I, ^=10.58 GeV, /37=0.425 


Case I, ^=3770 MeV 


Case II, ^/s=10.58 GeV, /37=0.425 


Decay Modes 


A v (xlQ-*) 


Av (xlO" 4 ) 


A v (xlO" 4 ) 


D+ K° s -k+ 


10.8 ( 9.0) 


15.9 (12.0) 


8.8 ( 8.5) 


D° -> Kl-n+K- 


12.9 (11.0) 


17.4 (14.7) 


10.5 (10.4) 


D° -» K%K+K- 


15.1 (12.8) 


30.6 (27.0) 


12.0 (11.8) 


B + _> K° s ir+ 


6.3 (4.5) 




5.2 (4.3) 


B° -> Kgjr+TT 


9.1 (7.1) 




7.5 (6.7) 


B° -> K° S K+K- 


9.5 (7.4) 




7.8 (7.0) 




2 4 6 -1 -0.5 0.5 1 
p K * (GeV/c) cose K * 



10"^ «v 

X 



a 
< 



10- 




0.5 



-1 -0.5 0.5 1 



p" 5 (GeV/c) 



cosG 



FIG. 3: Distributions of At> as a function of Kg momen- 
tum (left) and the polar angle (right) for D + — > K%-k + for 
Vs=10.58 GeV and /37=0.425 (top) and for D° -> K%-k + tt~ 
for -^5=3770 GeV configuration (bottom). Case I detector 
geometry is used in both instances. 



effect in very high energy experiments in the LHC envi- 
ronment can be smaller than the impact in experiments 
considered in this paper. 

In summary, we estimate the dilution effect in the mea- 
surement of CP asymmetry caused by the difference in 
nuclear interactions of K° and K° in e + e~ collisions for 
several typical experimental configurations. We find that 
the effect can be as large as 0.3% in decays involving 
low-momentum neutral kaons. The estimated systematic 
uncertainties on the calculated Ax> range in (15«30)% 
depending on K s momentum. We suggest that forth- 
coming high-sensitivity measurements of CP asymmetry 
involving neutral kaons in the final state should take into 
account the impact on the difference in K° and K strong 
interactions {A%>). 
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